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ABSTRACT 

Context. The stellar helium-to-metal enrichment ratio, AY/AZ, is a widely studied astrophysical quantity. However, its value is still 
not precisely constrained. 

Aims. This paper is focused on the study of the main sources of uncertainty which affect the AY/AZ ratio derived from the analysis of 
the low-main sequence (MS) stars in the solar neighborhood. 

Methods. The possibility to infer the value of the helium-to-metal enrichment ratio from the study of low-MS stars relies on the de- 
pendence of the stellar luminosity and effective temperature on the initial helium and metal abundances. The AY/AZ ratio is obtained 
by comparing the magnitude difference between the observed stars and a reference theoretical zero age main sequence (ZAMS) with 
the related theoretical magnitude differences computed from a new set of stellar models with up-to-date input physics and a fine grid 
of chemical compositions. A Monte Carlo approach has been used to evaluate the impact on the result of different sources of uncer- 
tainty, i.e. observational errors, evolutionary effects, systematic uncertainties of the models. As a check of the procedure, the method 
has been applied to a different data set, namely the low-MS of the Hyades. 
i-Q ■ Results. Once a set of ZAMS and atmosphere models have been chosen, we found that the inferred value of AY/AZ is sensitive to the 

ft' age of the stellar sample, even if we restricted the data set to low luminosity stars. The lack of an accurate age estimate of low mass 

field stars leads to an underestimate of the inferred AY/AZ of ~ 2 units. On the contrary the method firmly recovers the AY/AZ value 
for not evolved samples of stars such as the Hyades low-MS. Adopting a solar calibrated mixing-length parameter and the PHOENIX 

■+—* ■ GAIA v2.6.1 atmospheric models, we found AY/AZ = 5.3 ± 1.4 once the age correction has been applied. The Hyades sample 

—t ' provided a perfectly consistent value. 

Conclusions. We have demonstrated that the assumption that low-mass stars in the solar neighborhood can be considered as un- 
evolved, does not necessarily hold, and it may indeed lead to a bias in the inferred AY/AZ. The effect of the still poorly constrained 
efficiency of the superadiabatic convection and of different atmosphere models adopted to transform luminosities and effective tem- 

^ ' perature into colors and magnitudes have been discussed, too. 
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jy-! 1. Introduction actually measured helium abundance is not the original one, but 

fj ' the result of several complex processes, like dredge-up of nucle- 

g ■ It is a well known and firm result of stellar evolution studies that arly processed mate rial, diffusion and radiative levitation, which 

^ H ; the main structural, observational and evolutionary characters- severely alter the surface che mical composition. 

j. . tics of a star of given mass depend sensitively on the original As a consequencei in order to evaluate the original Y , the 

_ £ chemical composition, i.e. the initial helium and metal abun- only possibility is t0 rely on indirect met hods. This explains why 

ku dances, Y and Z, respectively. As a consequence, these parame- such an important param et eris still poorly constrained. 

£j ters affect also the observable quantities of stellar systems, from As egrly suggested by |Pei m bert& Torres-Peimbert I (fl97J 

03 star clusters to galaxies. from the analysis of the c h em ical composition of the HII regions 

While the present Z in a stellar atmosphere can be obtained in the Large Mag ellanic Cloud, a common approach in both stel- 

by the direct spectroscopic measurements of some tracer ele- lar and popu i ation synt hesis models is to assume a linear rela- 

ment, mainly iron, with the additional assumption on the mixture tionshin between the original Y and Z 
of heavy elements, the situation of Y is completely different. In 

the vast majority of stars the helium lines can not be observed, AY 

with the exception of those hotter than 20 000 K. This means that Y ~ Yp + ~/^ xZ ^ 
for the low-mass stars, which are the most common and long- 
lasting objects in the Universe, helium can be directly observed where Y P is the primordial helium content, i.e. the result of the 
only in advanced evolutionary phases, as the blue part of the hor- Big Bang nucleosynthesis, and AY/AZ the ratio which provides 
izontal branch or in the post-asymptotic giant branch. Thus, the the stellar nucleosynthesis enrichment. 

In the last four decades there has been a contin- 
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Indeed, as previously mentioned, the relationship between Y 
and Z adopted in stellar models directly affects some important 
quantities of stellar systems, both resolved and not, inferred by 
comparing observations and theoretical predictions. Thus, a pre- 
cise determination of lp and AY/AZ is of paramount importance 
for the studies not only of stellar evolution, but also of galaxy 
evolution. Furthermore, an accurate estimate of lp is of great 
cosmological interest, as it constrains the early evolution of the 
Universe, when the Big Bang nucleosynthesis occurred. In this 
paper we will focus on the value of the AY/AZ ratio. 

In the past, several techniques have been used to determine 
such a ratio by means of HII regions, both galactic and extra- 
galactic, planetary nebulae (PNe), the Sun and chemical evolu- 
tion models of the Galaxy (see Sect. [9] for a brief summary of 
these results). 

An alternative and well established way to determine the 
value of the helium-to-metallicity enrichment ratio takes ad- 
vantage of the dependence of the location of stars in the 
Hertzsprung-Russell (HR) diagram on their helium content. 
From the study of stellar populations in the galactic bulge 
Renzini ( 1994) inferred 2 < AY/AZ < 3. A frequently adopted 
approach relies on the analysis of the fine structure of the low- 
MS of the local field stars in the HR diagram. Pioneers of such 
an approach hav e been Faulkner, who found AY/AZ = 3.5 
dFaulknerlll96 7|) and Perrin and collaborators, who obtained 
AY /AZ = 5 (jP errin et al. 1977). Following these early stud- 
ies, |FernandeSjXebreton & Baglin ( 1996) constrained the value 
of AY/AZ to be larger than 2, by comparing the broadening in 
the HR diagram of the low-mass MS stars in the solar neigh- 
borhood and the theoretical ZAMS of several Y and Z. With 
a similar approach but taking advantage of Hipparcos data, 
IPagel & Portinari I (Il998l) inferred AY/AZ = 3 ± 2. This kind 
of app roac h culminated recently in the works bv lJimenez et al. I 
d2003l) and Casagrande et al. (2007) who provided, respectively, 
AY/AZ = 2.1 ± 0.4 and 2.1 + 0.9. 

The present analysis deals with the determination of the 
AY/AZ ratio by means of the comparison between the local K 
dwarf stars, for which accurate measures of both the [Fe/H] and 
parallaxes are available, and state-of-the-art stellar models. A 
great effort has been devoted to discuss the effect of the main 
uncertainties still present in stellar models on the inferred value 
of AY/AZ. 

In Sect. |2] we present the set of low-mass stellar models we 
have computed for this paper; in Sect. [3]we describe the data set 
we have used; Section|4]contains the description of the analysis 
method, while Sect. |5]deals with the possible sources of uncer- 
tainty that could affect the method itself. Results for the adopted 
data set are presented in Sect. [6] In Sect.|7]we investigate the pos- 
sibility of a non linear relation between Y and Z. In Sect. [8] we 
apply our method to an independent and unevolved set of stars, 
i.e. the Hyades low-main sequence. We compare our results to 



those obtained by other authors, with independent methods, in 
Sect. [9] SectionfTOlcontains the final discussion and summary of 
the whole paper. 



2. The models 

The stellar models have been computed on purpose for the 
present work with an updated version of the FRANEC evolu- 
tionary c ode which includes the s t ate-of the art input physics 
(see e.g. jC hieffi & Straniero 198 91: iDeglTn nocenti et al. 2008; 
Valle et al71 120091: iTognelli. Prada Moroni & DeglTnnocenti I 
2010). The main updating of the code with respect to previous 
versions include the 2006 release of the OPAL Equation of 



state (EOS) Q (see also Rogers et al. 1996) and, for tempera- 
tures higher than 12000 K, radiative opacity table s^! (see also 
Iglesias & Rogers 1996), while the opacities by Ferguson et al. 
(2005) are adopted for low er temperatures [f . The electron- 



conduction opacities are by Shternin & Yakovlev ( 2006) (see 
also Potekhin 1999). The opacity tables have been calculated 
by assuming the solar mixture by Asplund e t al.l (120051) . 

The extension of the convectively unstable regions is deter- 
mined by using the cl assical Schwarzschil d criterion. The mix- 
ing length formalism ( Bohm- Vitensel fl95 8 ) is used to model the 
super-adiabatic convection typical of the outer layers. As it is 
well known, within this simplified scheme, the efficiency of con- 
vective transport depends on a free parameter that has to be cali- 
brated using observations. We adopted the usual "solar" calibra- 
tion of the a of the mixing-length. More in detail, this means 
that we chose the value of a = 1 .97 provided by a standard so- 
lar model (SSM) computed with the same FRANEC code and 
the same input physics we used to compute all the other stellar 
tracks. 

Note that the "solar" calibrated value of the a param- 
eter strongly depends on the chosen outer boundary condi- 
tions needed to solve the differential equations describing the 
inner stellar structure, that is, the main physical quantities 
at the base of the photosphere (e.g. Montalban et al. 2004; 
Tognelli, Prada Moroni & DeglTnnocenti 2010). In order to get 
these quantities, we followed the procedure commonly adopted 
in stellar computations, which consists in a direct integration of 
the equations describing a ID atmosphere in hydrostatic equi- 
librium and in the diffusive approximation of radiative transport, 
plus a grey T(t) relationship between the temperature and the 
optical depth. Th e classical semiempirical T(r) relationship by 
Kris hna Swamv I ([T966) has been chosen. If a non-grey and more 
realistic model atmosphere is used, the solar calibrated value of 
a is different (see e.g. Tognelli, Prada Moroni & DeglTnnocenti 
2010, for a detailed discussion). 

Moreover, the "solar" calibrated value of the a parameter de- 
pends also on the input physics adopted in the SSM computation. 
Thus, to the sake of consistency, if the solar calibration approach 
is followed to fix the a parameter of a set of stellar models, the 
input physics and boundary conditions adopted to compute these 
models have to be the same as those used in the reference SSM. 

In spite of its widespread use, the solar calibration of the 
mixing-length does not rely on a firm theoretical ground, since 
there are not compelling reasons to guarantee that the efficiency 
of superadiabatic convective transport should be same in stars of 
different mass and/or in different stages of evolution (see e.g. 
Montalban et al. 2004, and references therein). However, for 
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what concerns the present paper, such an approach should be 
quite safe, since we deal with stars in the mass range 0.7-0.9 
Mq and that are on the Main Sequence. 

We also computed models with a different value of a, namely 
2.4, in order to quantify the effect of an uncertainty in the effi- 
ciency of the mixing-length on the inferred AY/AZ ratio. 

Our reference theoretical models have been transformed 
from the (log T e ff, log L/Lq) to the (B - V, My) diagram by 
means of synthetic photometry using the spectra database 
GAIA v2.6.1 cal culated from PHOENIX model atmospheres 
(Brott et aL 2005|). We per f ormed additional simulations adopt- 
ing the Castelli & Kurucz (2003) model atmospheres, to check 
the effect of the adopted color transformations on the inferred 
AY/AZ ratio. 

The original helium abundance in the stellar models has been 
obtained following Equation dT), where Z, once a solar mixture 
is assumed, is directly related to the observable [Fe/H] by 



6.0 - 



1-Fp 



1 T AZ 



(Z/X)q 



X 10-l Fe / H ] 



(2) 



We us ed fo r Y P the value 0.248, from llzotov. Thuan& Stasi 
(120071) and iPeimbert et aTl (120071) . a nd for (Z/XW the s olar 
metals-to-hydrogen ratio, 0.0165 from[Asplund et al. (2005). 

Stellar models have been calculated for 9 AY/AZ values (0.5, 
1, 2 ... 8) and 5 [Fe/H] values (from -0.6 to +0.2 in steps of 
0.2 dex). For each of these 45 combinations we calculated evo- 
lutionary tracks for 1 1 stellar masses (from 0.5 to 1 .0 Mq in 
steps of 0.05 Mq) in order to build zero age main sequence 
(ZAMS) curves that cover the whole HR region corresponding 
to the adopted data set (see Sect. 3). For each stellar mass we cal- 
culated its evolution starting from the pre-main sequence (PMS) 
phase. 

To determine the ZAMS position we used the following 
operative criterion: we calculated the local Kelvin-Helmoltz 
timescale for each model, i.e. ?kh = \Q\/L, where Q and L are 
the gravitational binding energy and the total luminosity; then 
we compared this number with the local evolutionary timescale, 

f ev = (j;^) , i.e. the inverse of the instantaneous rate of 
change of the luminosity. We found that a good operative def- 
inition of ZAMS is obtained by taking the first model for which 
r KH < 100 x f ev - 

As previously mentioned, an additional set of models, with 
the related ZAMS, has been computed for a value of the mixing- 
length parameter, a - 2.4. 

Thus, the present analysis can rely on a very fine grid of stel- 
lar models consisting of about a thousand evolutionary tracks 
calculated from the PMS phase to the central hydrogen exhaus- 
tion. 



3. The data set 

The stars of our samp le have been selected among the 
HIPPARCOS (lESA II 19971) stars with relative error on the paral- 
lax less than 5%. B and V band photometry are also taken from 
HIPPARCOS data set and they have typical errors between 0.01 
and 0.02 mag; combining them gives (B - V) colors with errors 
of the order of 0.03 mag. From the parallax and the observed V 
magnitude we computed the absolute magnitude My', with the 
quoted typical errors on V and parallax, the error on the abso- 
lute magnitude is ~ 0.1 mag, largely dominated by the error on 
the parallax. Given the small values of the distances, always less 
than 30 pc for the stars in our sample, we also assume that the 
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Fig. 1. The data color-magnitude diagram. Overplotted are three 
ZAMS all computed with AY/AZ = 2 and [Fe/H] values of -0.6 
(dotted), 0.0 (solid), 0.2 (dash-dotted). 



reddening is negligible, so that (B - V)o = (B - V). Our sample 
has been restricted only to those stars with My > 6 in order to 
take objects with long evolutionary time scales and hence mini- 
mize evolutionary effects. 

[Fe/H] determinations are taken from the Geneva- 
Copenhagen survey o f the Solar neighborhood catalog 
dNordstrom et all 12004 here after N04), f or 86 objects, and 
from the catalog by Taylor (Taylor 2005, hereafter T05) for 
21 objects. Of these 107 objects, 4 have [Fe/H] determinations 
from both catalogs; in these cases we only count the stars once, 
using their T05 metallicities. The total number is then reduced 
to 103 stars. N04 metallicity values are derived by a calibrated 
relation between Stromgren photometry measurements and 
spectroscopic determinations of metal abundances for a subset 
of objects. The T05 catalog is instead a collection of spectro- 
scopic determinations of metal abundances from the literature, 
where different results are put by the author on the same 
temperature scale; if more determinations are available for the 
same object, they are weight-averaged according to their quality. 
As pointed out in Tavlor (2005), determination of [Fe/H] from 
different authors may suffer strong systematic deviations, due to 
the different temperature scale chosen. The author showed that 
it is possible to reach a very good zero-point accuracy when 
data from different sources are put together in an appropriate 
way. Indeed a systematic difference in metallicity can be seen 
for the 4 stars that we have in common in our subsamples of the 
N04 and T05 catalogs. The metallicities for these objects are 
shown in Table Q] There we also show the new determinations 
of G eneva-Copenhagen met allicities using an improved calibra- 
tion dHolmberg et all 120071 hereafter H07). As shown by the 
authors, the N04 and H07 spectro-photometric calibrations give 
internally consistent results; only one of the 4 stars in common 
shows a small change in the [Fe/H] between the two catalogs. 
The average shift in [Fe/H] on this small sample and its standard 
deviation is<[Fe/H] T05 - [Fe/H] N04 >= -0.084 + 0.021 dex. 
Tavlor (2005) calcula tes the expected offse t between his 
metallicity scale and fhe lNordstrom et aTl (|2P04) one; according 
to his Table 10, for all the stars in our N04 sub-sample this offset 
is expected to be -0.023+0.017 dex, somewhat lower than what 
we find for our stars in common, which, anyway are only 4. We 
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Table 1. [Fe/H] for the 4 stars in common in the iTavlor I 
(2005) and Geneva-Copenhagen catalogs. T05, N0 4, H07 refer 
to ITav lor (2005), Nordstro m et all (120041) and Holmber g et al.l 
(2007) respectively. 



HIPPARCOS ID 


T05 


N04 


H07 


HIP 79190 


-0.410 


-0.340 


-0.300 


HIP 88972 


-0.157 


-0.080 


-0.080 


HIP 99711 


-0.006 


0.070 


0.070 


HIP 116745 


-0.336 


-0.220 


-0.220 



decided to apply the -0.023 offset to the Geneva-Copenhagen 
metallicities in their old version, i.e. the N04, since it is for the 
iNordstrom et all d2004l) calibration that lTavlorl d2005l) evaluates 
the zero-point offset. 

In the Geneva-Copenhagen survey the repeated radial veloc- 
ity measurements allow detections of almost all the possible bi- 
naries in the sa mple; we flagg ed out all the suspect binaries in 
the catalog. The iTavlorl (2005) [Fe/H] catalog is instead free of 
such contaminants, given its spectroscopic nature. 

Our final sample is made of 103 stars with good paral- 
laxes, photometry and [Fe/H] determinations. The typical errors 
are given by: cr(My) - 0.1 mag, <j{B - V) ^ 0.03 mag and 
cr([Fe/H]) si 0.1 dex. The error on the absolute magnitude is 
mainly due to the error on the parallax. 

Metallicities range from [Fe/H] = -0.6 dex to +0.2 dex; mag- 
nitudes range from My = 6.0 mag to My = 7.5 mag. The color- 
magnitude diagram (CMD) for the data is shown in Fig.Q] where 
the data are grouped in [Fe/H] bins 0.2 dex wide. 



4. The method 

The idea to determine the AY/AZ ratio using the position of low- 
mass MS stars in the HR diagram relies on the well known de- 
pendence of their luminosity and effective temperature on the 
original Y and Z. It is in fact firmly established that an increase 
of Y at fixed Z makes a star brighter and hotter. A decrease of Z 
at fixed Y leads to the same result. Such a behavior is the conse- 
quence of the effect on the opacity and mean molecular weight, 
i.e. the former gets higher as Z increases, while the latter grows 
with Y. 

As early shown by Faulkner (1967), who studied the effect 
of chemical composition variations on the position of theoreti- 
cal ZAMS, a simple but instructive explanation of this behav- 
ior can be obtained by means of homology relations (see also 
iFernandes. Lebreton & Baglin Ill9 96). Within this framework, it 
can be shown that varying Y and Z in such a way that AY » 5 AZ, 
leaves the bolometric magnitude M\, Q \ of ZAMS at fixed effec- 
tive temperature unchanged. This explains why the broadening 
of the local low-MS provides a AY/AZ indicator. 

However, this is not the whole story , as cl early proven 
by ICastellani. DeglTnnocenti & Marconi! dl999h . who com- 
puted a fine grid of full evolutionary models of ZAMS stars 
with several metal and helium abundances. They showed that 
the Mboi of ZAMS at a given T e ff depends quadratic ally 
on logZ and that only in a narrow range around the solar 
metallicity such a dep endence can be reasonably linearize d 
to AY/AZ- 5 (see also IFernandes. Lebreton & Baglinl Il996h . 
In addition, when comparing ZAMS models with real stars, 
one has to take into ac count the effect on the color indices. 
Castellani, DeglTnnocenti & Marconi (1 19991) showed that a 
AY/AZ w 7 is required to keep unchanged the B-V color at 
My-6 in a narrow range around Zq, The discrepancy between 



the AY/AZ required to keep the M^i and B-V unchanged clearly 
proves the crucial role played by the model atmospheres. In Sect. 
I5.4l we will further discuss the effect of different assumptions on 
the color-transformations on the inferred AY/AZ. 

The use of ZAMS models is allowed as long as observational 
data really lie on the ZAMS, or very close to it. Although this 
assumption is implicit in many studies which derive the AY/AZ 
from the fine structure of the low-MS, a detailed discussion on 
the effect of a deviation from such an assumption on the inferred 
value of the helium-to-metal enrichment ratio is still lacking. 
The reason is that very faint (My >6) local MS stars are usu- 
ally considered as if they were still on the ZAMS, since their 
evolutionary timescales are longer than the Galactic Disk age. 
We will show in Sect. l5.2l that this assumption is indeed critical, 
since underestimating the effects of age on both the position of 
the stars in the CMD and the diffusion of heavy elements below 
the photosphere leads to a severe bias in the final estimate of the 
enrichment ratio. We take this bias into account when giving the 
final result on AY/AZ for our work. 



Following Ijimenez et al. I (120031) . we did not directly use 
the broadening of the local low-MS in the CMD, instead, we 
compared models and data in a diagram like that of Fig. [2] 
right panel. After choosing a reference ZAMS, theoretical dif- 
ferences in magnitude, AMy, between that reference ZAMS and 
the other ZAMS curves, computed for different values of [Fe/H] 
and AY/AZ, are measured at a fixed value of the color index B-V 
(see Fig. |2l left panel). We checked that, within the current ac- 
curacy of the data, the derived AY/AZ value it is not affected by 
changing the reference ZAMS and/or the color index value. In 
fact, the ZAMS loci, in the range of magnitudes and colors that 
is involved here, are almost parallel to each other and the effect 
of the uncertain position of the star caused by the observational 
errors is much larger then that caused by a different choice of the 
reference ZAMS and/or the color index value. 

The differences AMy obviously depend on the chemical 
composition, i.e. on both AY/AZ and [Fe/H], as is clearly vis- 
ible in Fig. [2] (right panel). Observational differences between 
the data set and the same reference ZAMS are also measured in 
the way illustrated in Fig. [2] (left panel) and are plotted in Fig. [2] 
(right panel) as a function of [Fe/H]. 

To find the value of AY/AZ that gives the best fit to the data, 
we assign to each star errors in the three quantities My, B-V and 
[Fe/H]. The errors in My and B - V are assigned in the CMD, 
i.e. before the differences AMy are calculated; then the error 
in [Fe/H] is assigned in the (AMy, [Fe/H]) diagram. Magnitude 
and color errors are considered to be distributed as gaussian with 
<x equal to the quoted uncertainty for that star; the assumption of 
gaussian errors is reasonable, considering that they come from 
the HIPPARCOS photometric errors plus (for the absolute mag- 
nitude) the HIPPARCOS parallax error; these two sources of er- 
ror are independent and our objects are all close by and have well 
determined parallaxes, so that they don't suffer the Lutz-Kelker 
bias. Regarding the [Fe/H] value the situation is different, since 
the errors associated to each value are highly affected by system- 
atic effects like the choice of the temperature scale; for this rea- 
son, and since we can not reconstruct the real error distribution 
in [Fe/H], we adopted an uniform distribution of [Fe/H] errors. 
Anyway, we have also checked that a gaussian distribution for 
[Fe/H] values leaves the results essentially unaffected. 
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Fig. 2. Illustrative example of how theoretichal and observational differences are computed and used to produce a AMy vs. [Fe/H] 
diagram. Left: Theoretical differences are obtained with respect to a reference ZAMS computed using AF/AZ = 2 and [Fe/H] = 0.0 
dex (solid line) at the reference color of B - V - 0.95 (dotted vertical line); the dashed lines are ZAMS computed for AF/AZ = 4 
and [Fe/H] e [-0.6, 0.2]. Observational differences (represented by the distance between filled and empty stars) are calculated for 
an observed star at a given color with respect to the same reference ZAMS. Right: The differences are used to place both ZAMS 
and stars in the AMy vs. [Fe/H] diagram. The solid line corresponds to the curve of AF/AZ = 4 (the ZAMS of the left panel); the 
dot-dashed line to AF/AZ = 0.5; the dotted line to AF/AZ = 8. 



Once a "new" data set is created from the original value plus 
the errors, we determine the theoretical curve which minimizes 
the quantity: 



*%-E 



AM VJ - AM K/ ([Fe/H],) 



o-(AM V j) 



(3) 



where j runs over the curves (i.e. over different AF/AZ values) 
and i over the data. The quantity AMy./fFe/H],-) is the value of 
the magnitude difference for the j-th theoretichal ZAMS, calcu- 
lated for the [Fe/H] value of the j-th star. The scheme is iterated 
in order to obtain not only the best fitting value but also an esti- 
mate of the error on the final result due to the observational un- 
certainties. For each iteration we take the value that minimizes 
Equation ([3]) and then plot an histogram of the number of best-fit 
occurrences for each AF/AZ in our models grid. A gaussian fit 
to the histogram is performed and the mean and the standard de- 
viation of this gaussian are used as estimators of the true AF/AZ 
and its error. Figure [3] shows the histogram of occurrences for a 
Monte Carlo simulation with a total of 10 5 iterations using our 
data set of 103 stars; superimposed is the best-fitting gaussian. 
We will discuss in Sect. [6] why, in this particular case, a gaussian 
curve does not fit well the histogram of occurrences, with the 
central part of the histogram which is too broad and flat to be 
well approximated by a single gaussian. We just anticipate that 
this discrepancy is mainly due to the offsets in the [Fe/H ] scales 
for N0 4 and T05, even after the correction of Table 10 of lTavlorl 
(2005) is applied to the N04 data. 



5. Analysis of possible sources of uncertainty 
using artificial data sets 

To check the reliability of the procedure described above and 
to evaluate the contribution of the various possible sources of 
uncertainty, we applied our method to a number of artificial data 
sets with controlled input parameters. 
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Fig. 3. Results of 10 5 simulation runs to determine the best fitting 
AF/AZ value for our 103 stars data set, according to Equation 
0. The best-fitting gaussian is also shown, together with its 
mean value and standard deviation. 



These sets have been created by interpolation in our fine grid 
of stellar models. Stellar masses are randomly generated from a 
power law initi al mass function (IM F) ^ = c x nT a with a slope 
of a - 2.3 (ISalpeter Hl955tlKroupal20oi1) : an IMF with a unique 
value of the exponent is a good law in our range of simulated 
masses (0.5 < M s [ m /MQ < 1). The chemical composition is 
calculated by fixing the input value of (AF/AZ) in and extracting 
random values for [Fe/H] for each star; F and Z are then calcu- 
lated using the two Equations ([TJ and (0. In what concerns the 
stellar ages generation, we performed different kind of simula- 
tions, adopting three age-laws (i.e. star formation rates, SFRs), 
namely a Dirac's delta centered about a given age (i.e. coeval 
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stars), a uniform and an exponentially decaying age distribution, 
respectively, in the range 0-7 Gyr, a reasonable approximation 
for the age of the galactic disk. 

Since the FRANEC code includes a treatment of diffusion 
of the elements, given a star of age t, we take from the models 
the corresponding surface value of [Fe/H] T , which is different 
from the initial value [Fe/H]o because of diffusion itself. This 
is the value that we use for the recovery, because in a real star 
the observed [Fe/H] is the present one and not the initial. Note 
that neglecting diffusion may lead to a bias (an underestimate, 
indeed) in the final AY/AZ if the sample of stars is old enough to 
have experienced a not negligible amount of diffusion of heavy 
elements; we will show this in Sect. 15.21 

When generating the stellar models parameters for our sets 
of stars, we do not take into account any age-metallicity relation, 
hence [Fe/H] and age values are independ ently and ra ndomly ex- 
tracte d. Recent works ( Nordstro m et al J [20 04; Holmberg et alJ 
2007) have indeed shown that there is no evidence of an age- 
metallicity relation for local disk stars. Given mass, age and 
chemical composition, we interpolated in our fine grid of stellar 
models to obtain the observational properties of the simulated 
stars. The number of stars in each simulated set is 1 10, a number 
comparable to the 103 stars of our real data sample. 

In the following and in the related figures, we will refer to the 
input parameters of our simulated samples using the subscript in 
and to the output of the recovery method using the subscript out . 

5.1. The effect of measurements errors 

The first test we have performed was made to check whether our 
recovery method was able to get the right AY/AZ from an "ideal" 
sample of stars affected only by observational errors on the mag- 
nitude, color and [Fe/H] values. By ideal we mean a sample that, 
regardless of masses and [Fe/H] distribution (which indeed were 
generated in a completely random fashion), contains only stars 
really lying on the ZAMS. It is worth to point out that in the case 
of the real data, this is only a simplifying assumption, which 
can not be exactly fulfilled, since the observed stars in our sam- 
ple have unknown ages that span the whole range of ages in the 
Galactic disk. 

Once the artificial sample has been generated, we associated 
to each star an error in absolute magnitude, color and [Fe/H] 
typical of our real sample of data, i.e., 0.1 mag, 0.03 mag and 
0. 1 dex respectively (see Sect. [3} ■ 

We found that our recovery method is not affected by obser- 
vational errors of this order of magnitude. As it is possible to 
see in Fig. [4] given a (AY/ AZ) in of 4, the best value that comes 
out from our Monte Carlo method and the gaussian fit to the his- 
togram of occurrences is indeed (AF/AZ) out = 3.87 + 0.63. So 
the outcome of the method is perfectly consistent with the input 
value and; moreover it has a very small range of 1 .26 at a level of 
lcr, which is comparable to the resolution in our models grid (1 
unit), and which we may quote as the nominal or intrinsic error 
of the method, associated to the typical error of the actual data. 



5.2. The effects of age and heavy elements diffusion 

Although our stellar sample has been obtained by selecting very 
faint stars (My > 6 mag, i.e. M < Q.9Mq, the actual value de- 
pending on the chemical composition), we found that evolution- 
ary effects strongly affect the final result introducing a non neg- 
ligible bias. This is actually one of the most important results of 
this work. Thus, one should be very careful in properly taking 
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into account the evolutionary effects, i.e. the displacement from 
the ZAMS, when helium- to-metals enrichment ratio is derived 
from the low MS fine structure. 

Indeed, by creating artificial stellar data sets with 
(AF/AZ); n = 4, but no longer on their ZAMS position, we found 
an output value of our recovery method of (Ay/AZ) out ~ 2. 
The actual best-fit value for each simulated data set depends 
on the exact parameters used to generate the artificial sample, 
like the maximum age or the functional form of the age distribu- 
tion (uniform or with an exponentially decaying SFR). The total 
number of stars, their [Fe/H] and magnitude ranges are always 
kept the same between the different simulated data set. Figure 
|5] shows the results of our method using data sets where the 
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Fig. 6. The evolution of stars mimics lower values of the en- 
richment ratio. A 7 Gyr isochrone calculated for AY/AZ = 4 
and [Fe/H]=0 is very similar to a ZAMS with AY/AZ = 2 and 
[Fe/H]=0. The arrow indicates the effect of evolution on the po- 
sition of the isochrones. 



Fig. 7. Effect of changing the mixing length parameter from a = 
1.97 to 2.4 on the calculated ZAMS; all the models shown are 
calculated for [Fe/H] = +0.2 dex. The arrows indicate in which 
direction the ZAMS move in the CMD when increasing a (upper 
arrow) or AY/AZ (lower arrow). 



simulated stars have ages uniformly distributed between and 
7 Gyr. The dashed line indicates the results when the [Fe/H] val- 
ues associated to each star at a given age are the same as the 
ones at the ZAMS; whereas the solid line indicates the results 
when diffusion is taken into account, i.e. the value at a given 
time [Fe/H], is different from [Fe/H]zAM5- Even without tak- 
ing into account diffusion, the method gives an output value of 
(AY/AZ) our = 2.24 + 0.51, quite different from (AY/AZ) in = 4. 
When also diffusion is taken into account, the best fit value suf- 
fers an additional shift, with (AY/AZ) out - 1.81 ± 0.63. Since 
we don't know what is the real age distribution of the stars of 
our sample we can not really quantify the bias, but after many 
experiments with several data set, we conclude that it must be of 
the order of S age (AY/AZ) = -2 ± 0.5. 

The effect of evolution on the derived enrichment ratio is 
easy to understand by looking at Fig. [6] Here two ZAMS 
with different values of AY/AZ, namely 2 and 4, and same 
[Fe/H] = 0.0 are shown together with isochrones of 1 and 7 Gyr, 
calculated with AY/AZ - 4 and [Fe/H] = 0.0; it is clear that evo- 
lution causes a shift of the whole curve towards redder colore in 
a completely indistinguishable fashion as a lower AY/AZ does. 

As a summary, this means that, even if our observational data 
set has been selected with a very strict cutoff of My - 6 mag, 
evolutionary effects still play an important role. The real, unbi- 
ased value of AY/AZ coming out from our analisys has then to 
be corrected, by subtracting $ age from the nominal value given 
by the Monte Carlo method. We then expect that the real value 
is higher by about two units than what can be found by blindly 
applying this method to the data. 



5.3. The effects of the uncertainty on the mixing-length 
efficiency 

The current generation of stellar models is not yet able to firmly 
predict the effective temperature of stars with a convective en- 
velope, such as those belonging to our sample. The reason is 
that a satisfactory and fully consistent theory of convection in 



4 A single star becomes hotter and more luminous after leaving the 
ZAMS, here we are referring to the overall shape of the isochrones and 
ZAMS in this range of magnitudes and colors. 



superadiabatic regimes is still lacking, hence a very simplified 
approach is usually followed. The approach commonly adopted 
in the vast majority o f evolutionary codes is to implement the 
mixing-length theory (Bohm-Vitense 1958), in which the aver- 
age efficiency of convective energy transport depends on a free 
parameter a that must be calibrated. Our reference set of stellar 
models has been computed adopting our solar calibrated value 
of the mixing length parameter, namely a — 1.97. Nevertheless, 
we calculated a whole new grid of models using a = 2.4 in or- 
der to evaluate the effect of this still uncertain parameter on the 
derived value of AY/AZ. As well known, the predicted effective 
temperature of a stellar model with a convective envelope is an 
increasing function of the value of the mixing length parameter 
a, as a consequence of the shallower temperature gradient due 
to a more efficient convective energy transfer. 

Figure|7] shows the effect of different adopted a values on the 
calculated ZAMS. The assumed mixing-length parameter affects 
also the inferred AY/AZ ratio, since it directly influences the pre- 
dicted position in the HR diagram of the ZAMS models. As one 
can easily see in Fig. [7] in order to recover the ZAMS locus of 
models computed with an higher value of a, a lower AY/AZ ra- 
tio is needed. Notice also that, the impact of the mixing-length 
efficiency on the predicted effective temperature of ZAMS mod- 
els of the same chemical composition becomes progressively 
smaller at faint magnitudes, i.e. for very low-masses. Such a be- 
havior is the consequence of the almost adiabatic nature of con- 
vection in the envelopes of very-low mass stars (M < Q.IMq), 
characterized by high densities and low temperatures. 

Since an increase in a affects the models in the same direc- 
tion as an increase in AY/AZ, we find, as expected, that our re- 
covery method gives a lower value of the enrichment ratio when 
ZAMS calculated for a - 2.4 are used, that is, (AY/AZ) out - 
2.00 + 0.61, quite different from (AF/AZ),„ = 4. This is shown in 
Fig. [8]which reports the output of the Monte Carlo method, when 
the same data set of Sect. 15. H is used, i.e. a set of stars lying on 
the ZAMS calculated with a = 1.97. We already mentioned that 
the use of the solar calibration should be the safer choice when 
dealing with Main Sequence stars of solar-like masses such those 
of our data set. On the other hand, this numeric experiment al- 
lows to quantify the effect of a wrong assumption of the mixing 
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our Monte Carlo method using ZAMS calculated with a mixing 
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Fig. 9. ZAMS transformed into the observational plane using 
PHOENIX (solid lines) and ATLAS9 (dashed lines) model at- 
mospheres. All the ZAMS shown have been calculated with 
AY/AZ = 4 and [Fe/H] = -0.6 (left) or +0.2 (right). 



Glass photometric system ze ro points are those of Table Al in 
iBessel. Castelli & Plez I (Il999l) . 

Both model atmosphere grids completely cover our range of 
T e g, logg and [Fe/H]. An exhaustive description of these models 
is far beyond the scope of this paper, nevertheless we want to 
mention briefly some of the most important differences between 
them. The ATLAS9 code solves the radiative transfer equation 
in a plane parallel atmosphere, while PHOENIX code takes into 
account the curvature of the atmosphere, even though in spheri- 
cal symmetry, i.e. keeping a ID approach. This difference is not 
very important for dwarf stars, for which the curvature (i.e. the 
ratio of the extent of the atmosphere to the radius of the base 
of the atmosphere) is very small. More important is the differ- 
ent database of molecular opacities for the two models. While 
atomic opacities databases are similar between the two models, 
PHOENIX models include a lot of molecular species (~ 650) 
and molecular transitions (~ 550 millions) which play an impor- 
tant role for low mass stars, specially as the metallicity increases. 

Figure [9] shows the comparison between ZAMS models 
transformed using the two sets of model atmosphere. It is evi- 
dent how the PHOENIX ZAMS are always redder in the B - V 
color than the ATLAS9 ones at fixed magnitude. Moreover the 
differences in stellar colors increase with decreasing mass (in- 
creasing magnitude) and also, at a given magnitude, the differ- 
ences increase with increasing metallicity (from left to right in 
Fig. [9]). Note that both a lower temperature and a higher metal 
content favor the formation of chemical composites; in particular 
the first molecules start forming when the effective temperatures 
drops below ~ AQQQK. 

The sizeable difference in magnitude and color index be- 
tween the same theoretical ZAMS transformed into the obser- 
vational plane by the two quoted model atmospheres, directly 
translates into a large difference in the inferred AY/AZ value. 
We used the ATLAS9 ZAMS, running our recovery method on 
a simulated data set with (AT/AZ) in = 4 but generated using the 
PHOENIX ZAMS. The ATLAS9 ZAMS are so much bluer than 
the PHOENIX ones that in each iteration of the method we al- 
ways find the lowest possible value of (AF/AZ) out available in 
our grid of models, i.e. 0.5. 

The uncertainty due to the chosen model atmosphere is then 
by far the most severe source of uncertainty affecting the final 
value of the enrichment ratio, at least among the uncertainties 
coming from the models side. 



length parameter in the adopted stellar models on the inferred 
helium-to-metals enrichment ratio. 



5.4. The effects of different transformations from the 
theoretical to the observational plane 

To compare observational data with stellar models, one needs to 
transform the theoretical predictions of evolutionary codes from 
the (log Tgg, log L/Lq) plane (HR diagram) to the observational 
plane, in our case the (B - V, My) Color-Magnitude diagram. A 
common procedure is to use synthetic stellar spectra, calculated 
using model atmosphere codes, and to convolve them with the 
filter throughput of the photometric system needed; this tech- 
nique, referred to as synth etic photometry is described in detail 
in, e.g. Girardi et al. (2002). We transformed our models using 
PHOENrx TBrott et alj|2005h and ATLAS9 dCastelli & Kurucz I 
2003) model atmospheres; the UBVRIJHKL Johnson-Cousin- 



5.5. The effects of different choices for the heavy elements 
mixture 

As previously explained, all the models in our grid have 
been c alculated using the solar-scaled mixture bv lAsplund et al.l 
(2005), hereafter AGS05. A new version of the solar mixture 
has be en recently published by the same group in Asplund et al. 
(I2009h . hereafter AGSS09. The true dependence of the inferred 
AY/AZ on the mixture choice could be evaluated only by re- 
calculating an equivalent grid of models using the new mixture 
and by re-running the whole procedure illustrated in this work. 
However, the exact determination of the solar mixture is still 
an open problem, thus a recalculation of all the models is not 
needed, in our opinion, until this issue will be definitively set- 
tled. Nevertheless, to have an idea of the influence of a variation 
of the solar mixture on our results, we computed two new sets 
of ZAMS with different mixtures and compared them with our 
reference AGS05 ZAMS in the CMD. This gives at least an in- 
dication of how the inferred enrichment ratio may depend on the 
mixture. In addition to the ZAMS for AGS05 and AGSS09 mix- 
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Fig. 10. ZAMS calculated for the three different mixtures, 
GN93, AGS05 and AGSS09, with AY/AZ = 4 and [Fe/H] = - 
0.6 (lower B-V values) or [Fe/H] = +0.2 (higher B-V values). 



tures, we also com puted and compared ZAM S calculated with 
the older mixture bv lGrevesse & Noels] d 1993b . hereafter GN93, 
still widely used in the literature. We highlight the fact that 
we can control the effect of the mixture changes on the stellar 
structure, by using opacity tables calculated with different mix- 
tures both for the high- temperature (OPAL) and low-temperature 
dFerguson et al.l d2005l) ) opacities. On the other hand we can not 
evaluate the effects of different mixtures on the model atmo- 
sphere, i.e. on the transformation of our theoretical models from 
the theoretical HR diagram to the CMD, since PHOENIX syn- 
thetic spectra are available only for a single mixture of heavy 
elements. 

The effect of the adopted solar mixture on the models is 
twofold. First, for a given global metallicity Z, changing the in- 
ternal distribution of metals will mainly affect the opacity and 
the nuclear burning efficiency (via the CNO abundances). For 
models with the same Y and Z, the ZAMS computed with the 
GN93 mixture are bluer than that with the AGSS09 one, which, 
in turn, is bluer than the ZAMS computed with the AGS05 
mixture (see e.g. DedTnnocenti, Prada Moroni & Ricci 2006; 
Tognelli, Prada Moroni & Degl'Innocenti 2010). Second, for a 
given choice of the two values of AY/AZ and [Fe/H], it is clear 
from the coupled equations (Q]) and (f2]i that the values of Y and 
Z will be different, given the different values of (Z/X)q. More in 
detail, the GN93 mixture provides models with higher Z. These 
models are thus redder on the CMD, than those calculated with 
AGSS09 mixture. The latter mixture, in turn, provides models 
that have higher Z and appear redder than the AGS05 mixture. 
Thus, the effects on the opacity and on the scaling, i.e. the ac- 
tual value of Z given [Fe/H] and the different (Z/X)q, affect the 
ZAMS position in opposite directions. 

We calculated models for the two additional mixtures by fix- 
ing AY/AZ = 4 and choosing the two values [Fe/H] = -0.6 and 
0.2. The results of the calculations are shown in Fig. \W\ The 
corresponding Y and Z values for each mixture are shown in 
Table|2] From Fig.[10]one can see that GN93 models seem to be 
slightly redder than our reference AGS05 models; when com- 
pared to observational data, an higher value of AY/AZ would 
then be needed to reproduce observations, as compared to the 
AS05 mixture. AGSS09 ZAMS go in the opposite direction 
meaning that the inferred value of the enrichment ratio would 
be lower when compared to AGS05. In the comparison between 
AGS05 and AGSS09, which have quite similar (Z/X)q, the ef- 



Table 2. Y and Z values for the three different adopted mixtures 
(see text in Sect. 15.5b . for two reference values of [Fe/H] and 
AY/AZ - 4. The solar (Z/X)q values corresponding to the three 
mixtures are also indicated. 





GN93 

(Do = - 0245 


AGS05 

(Do = o- 0165 


AGSS09 

(Do = - 0181 


[Fe/H] = 


-0.6 


Y = 0.2660 
Z = 0.00449 


Y = 0.2602 
Z = 0.00305 


Y = 0.2614 
Z = 0.00334 


[Fe/H] = 


0.2 


Y = 0.3458 
Z = 0.0245 


Y = 0.3178 
Z = 0.0174 


Y = 0.3235 
Z = 0.0189 



feet on the opacity and burning efficiency seems to prevail on the 
scaling; viceversa in the comparison between AGS05 and GN93, 
the huge difference in (Z/X)q compensates for the opacity effect, 
and the much higher Z value in the GN93 case brings the ZAMS 
models toward redder colors when compared to AGS05. 

The relative behavior of different ZAMS is moreover a non 
trivial function of the actual value of [Fe/H], given the non linear 
relation between the two couples of values (AY/AZ, [Fe/H]) h> 
(Y, Z), represented by the two equations (Q]) and (f2j). Overall it 
seems that changing the mixture affects our method in a non 
negligible way, but a totally consistent check could be done 
only when the appropriate model atmospheres are used to do 
the transformation from the HR diagram to the CMD and only 
calculating complete grids of models for different mixtures. 

6. Results using the observational data 

After having carefully checked the capability of our recovery 
method and having studied many uncertainty sources by means 
of controlled artificial data sets, we applied the above technique 
on the sample of real observational data described in Sect. [3] 

Figure [3] shows the result of the Monte Carlo recovery 
method applied to the local low-MS field stars, which provides a 
nominal value of the enrichment ratio of AY/AZ = 3.32 + 1.29. 

This result for the enrichment ratio is obtained using our 
favorite set of models, i.e. ZAMS calculated with the solar 
a - 1.97 and transformed using the PHOENIX model atmo- 
spheres. We showed in Sects. 5.3 and 5.4 the significative effects 
on the inferred helium-to-metals enrichment ratio of a different 
choice of the mixing length parameter a and of the model atmo- 
sphere, respectively, in the case of synthetic data. The recovery 
procedure applied to the real observational data provides a nomi- 
nal value of the enrichment ratio of AY/AZ = 1.59 + 1.01 adopt- 
ing the set of theoretical models computed with a - 2.4 and 
transformed into the observational plane by means of PHOENIX 
atmosphere models and AY/AZ - 0.5 adopting our standard set 
of models with a - 1 .97 but the ATLAS9 atmosphere models. 

Moreover, we will show a self-consistency check in Sect. [8] 
to prove that inside our choice the result we obtain for AY/AZ 
can be used to calculate models that fit very well an independent 
data set, in particular the Hyades main sequence. 

In Sect. [3] we mentioned the fact that the fit of a gaus- 
sian distribution to the histogram of occurrences of Fig. [3] is 
not perfect. The actual distribution looks quite flat in the range 
AY/AZ e [2, 4] and the <x of the fitting gaussian is twice as 
large than what we would expect from the observational un- 
certainties alone (see Sect. 15. It . We also mentioned that our 
[Fe/H] estimates come from two different catalogs, and that the 
N04 data have been re-zeroed according to T05 transformations. 
Nevertheless, when considering the two subsets of data sepa- 
rately, a large difference in the final result is still visible, as if the 
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M0 H 



(A Y/A Z) oul = 2.96 ± 1.44 (T05) 
(A Y/A Z) m = 4.51 ± 0.56 (N04) 




Fig. 11. Results of our Monte Carlo method when applied sep- 
arately to the subsets of objects with [Fe/H] from T05 and N04 
(the latter corrected using Table 10 of T05). The solid lines are 
for T05, the dashed for N04. Overplotted are the best-fit gaus- 
sians. 

different [Fe/H] scale adopted, even after the re-zeroing, would 
mimic a bimodal chemical composition distribution. From Fig. 
[TTl it is clear that the two results for the separate subsets are in- 
consistent at lcr level. The final result we get for the whole set, is 
then an average of the two distributions of Fig.QT] unfortunately 
this method is quite sensitive to the adopted scale of [Fe/H] de- 
terminations. We won't choose any of the two subsets as the best 
one, we just state that systematics is probably still affecting the 
[Fe/H] of our objects, even after the correction of N04 zero point. 
Anyway, in order to proceed further with the analysis we will 
use the results of the gaussian fit to the histogram of occurrences 
for the whole sample of data as our nominal value for AY/AZ. 
As we showed in Sect. 15.21 to obtain the real AY/AZ value, the 
nominal value of AY/AZ - 3.32 + 1.29 provided by the recov- 
ery method must be corrected for the age-bias, i.e. the artificial 
shift in the recovered value of the enrichment ratio caused by 
neglecting evolutionary effects and diffusion. We estimated the 
extent of such a correction to be of S a 



ge 



-2 + 0.5 by means of 
Monte Carlo simulations. Thus, the final corrected value of the 
helium-to-metals enrichment, combining the errors in quadra- 
ture is AY/ AZ = 5.3 + 1.4. 

7. A possible non-linear relation between Y and Z? 



In their recent work. lCasaprande et alj (2007) suggested that the 
assumption of a linear relation between metals and helium en- 
richment of the ISM, which leads to the definition of the enrich- 
ment ratio itself (see Equation Q]) may not necessarily be valid. 
We tried to check whether it is possible to unambiguously in- 
fer such a non linear behavior from the presently available data 
accurateness. We performed the test by means of an artificial 
and controlled data set built interpolating our stellar models. The 
synthetic stellar sample mimics the main characteristics of the 
real data set. In particular, we simulated stars divided in metal- 
licity bins 0.1 dex wide in the range [Fe/H] e [-0.6, +0.2], with 
the same number of stars in each bin as our observational data, 
for a total number of 103 objects. To avoid the additional diffi- 
culties introduced by the age-bias, the simulated stars were taken 



on their ZAMS position. We then ran our Monte Carlo method 
bin by bin to see whether our recovery is able to recover the in- 
put value for AY/AZ in each bin, given the typical observational 
errors on our objects. 

FigureQ~2]shows the outcome of such an experiment. On the 
two panels each point represents the result of the Monte Carlo 
method when applied to the corresponding [Fe/H] bin, with its 
lcr uncertainty; both values come from the gaussian fitting of 
the histogram of occurrences. The dashed lines represent the re- 
sults of the gaussian fit the histogram of occurrences when the 
method is applied to the whole sample of simulated stars; the 
shaded area corresponds to the lcr uncertainty around the mean. 
The two panels represent results for two different ways of bin- 
ning the data; on the left panel we show the results when data 
are binned in 0.1 dex wide non overlapping bins; on the right 
panel the results are for 0.2 dex wide bins, in this case adjacent 
bins have an overlap of 0.1 dex, meaning e.g. that a [Fe/H] = 
-0.35 star can be found in both the [-0.5, -0.3] and the [-0.4, - 
0.2] bins. Regardless from the way the binning is performed, it 
follows from this numerical experiment that the method works 
fine only for [Fe/H] values equal or greater than -0.3 dex. AY/AZ 
determinations with this method, considering only lower values 
of the metallicity are unreliable. The missing points in the two 
panel, correspond to bins where the gaussian fit failed, simply 
because the histogram of occurrences is too irregular to be fitted 
by a gaussian. 

Hence any discussion on a possible non linearity of the 
helium-to-metals enrichment relation, extending over a wide 
range of metallicities, can not be carried on when Y and Z value 
are determined from low-main sequence fitting. The fine struc- 
ture of the low main sequence at values of the metallicity lower 
than -0.3 dex is indeed too weakly dependent on the chemical 
composition, to be used to infer any non linear behavior of the 
ISM enrichment in helium and metals. This is quite clear from 
the right panel of Fig. [2] which shows that the separation among 
curves with different AY/AZ becomes progressively smaller and 
smaller as the [Fe/H] decreases. At low enough values of [Fe/H], 
this separation is simply too small with respect to the spread of 
the observed point, due to the current uncertainties in both mag- 
nitudes and [Fe/H] , which are of the order of 0. 1 mag and 0. 1 dex 
respectively, to allow a correct determination of the stars chem- 
ical composition. This uncertainty is reflected in the large error 
bars of the lower metallicity bins of Fig. [12] 

We are not claiming that any non linear law for the helium- 
to-metals enrichment of the ISM must be ruled out; we are just 
stating that better data are needed to confirm this possibility by 
making use of the fine structure of the low main sequence. The 
GAIA mission will provide us with better parallaxes, which will 
yield better absolute magnitude determinations. Advances in 3D 
and NLTE stellar atmospheres modeling will also produce bet- 
ter constrained spectroscopic abundances determinations in the 
future. 



8. The Hyades test 

As previously discussed, the AY/AZ provided by the recovery 
procedure must be corrected for the age-bias, since the above 
described technique uses the theoretical ZAMS loci as reference 
mile stones in the CMD. Hence, any unaccounted displacements 
of the observed stars from the ZAMS due to evolutionary effects 
translates directly in an underestimate of the inferred AY/AZ. 
Clearly, the younger the stellar sample analyzed, the smaller 
such unaccounted displacements and the smaller the related age- 
bias. 
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Fig. 12. Output values of the recovery method when data are binned in metallicity (see text for details). The green dashed lines 
represent the (AY/AZ) 0Ut when the whole simulated sample is used. The input value is (AF/AZ); n = 4. 



To safely neglect any possible bias due to evolution, 
we also applied our recovery method to a very young stel- 

> 5.2) of the 
age of about 



lar sample: the stars of the lowest part (M\ 
Hyades MS. For this well studied cluster an 



500 f 600 Myr is generally estimated (see e .g. Perryman et 



1 998t ICastellani, Peel' Innocenti &P rada Moroni I 2001 



in 
m 



Castellani et al. I 120021) . Thus, the selected low luminosity 



stars, which correspond to masses lower than about 0.9Mq, 
are really on, or at least very close to, their ZAMS position. 
Moreover, the data available for this cluster are even more 
precise than those of our data set for field stars in the solar 
neighborho od; using the kinematically correcte d parallaxes 
data from iMadsen. Pravins & Lindegren] (120021) the typical 
errors are cr(My) - 0.05 mag and cr(B - V) =* 0.02 mag. 
The metallicity value we use for the Hyades is taken from 
IPerrvman et all (Il998l) : [Fe/H] = 0.14 ± 0.05 dex. 

By running our Monte Carlo method, the value found for 
the enrichment ratio is {AY / AZ)n\ad - 4.75 + 0.35, in perfect 
agreement with what we found using local MS stars and taking 
into account the evolutionary bias. 

To further check the consistency of this final result, we cal- 
culated isochrones using the closest value of the enrichment ra- 
tio in our grid, i. e. AY /AZ = 5, and for the extreme values of 
IPerrvman etaD (fl998) [Fe/H] interval, i.e. 0.09 and 0.19 dex. 
Models are calculated for our solar calibrated value of the mix- 
ing length parameter and transformed to the observational plane 
using PHOENIX model atmosphere. Average extinction towards 
the Hyades is negligible and this cluster is old enough not to 
present any intra-cluster material, so also differential extinction 
can be ignored; hence we assume (B - V)o — (B - V) b s . 

Figure [T3l shows the CMP of the Hyades with overimposed 
our isochrones. The good agreement between our models with 
the chemical composition derived by the recovery method and 
the Hyades data is quite encouraging and make us confident 
about the adopted procedure. Even more important is that the 
fit to the Hyades MS is good also at higher luminosities. Stars 
with M > \.2Mq (corresponding to a magnitude My < 3.6, the 
exact value depending on chemical composition and age) have 
essentially radiative envelopes, meaning that their predicted tem- 
peratures are not affected by the choice on the mixing length 
parameter, a. 



(Y, Z) - (0.3116, 0.0136) 

(Y, Z) = (0.3313, 0.0167) 




Fig. 13. The CMP for the Hyades Main Sequence. Isochrones 
for two values of the age are overplotted: 500 Myr (black lines) 
and 600 Myr (red lines). Pashed and solid lines correspond to 
two different values of the chemical composition. Y and Z values 
are obtained from Equations ([TJ and (O using [Fe/H] = 0.09 dex 
(dashed lines) and [Fe/H] = 0.19 dex (solid lines). 



9. Comparison with independent methods 

Many independent approaches have been followed in the past 
to determine the helium-to-metals enrichment ratio. A detailed 
comparison between the value of helium-to-metal enrichment 
ratio AY/AZ obtained by means of different techniques is beyond 
the aim of the paper, since it would require a careful discussion 
of the different uncertainty sources, both systematics and ran- 
dom, affecting the various techniques. Furthermore, the primor- 
dial helium abundance Y p adopted by different authors and/or 
in different periods can be very different, affecting the inferred 
AY/AZ value. On the other hand, we think that a brief descrip- 
tion of the results provided by different approaches might be of 
interest. 

One of the most fruitful consists in observing HII re- 
gions in the Milky Wa y and in other galaxies. The p ioneer- 
ing papers by[P eimbert & To rres-Peimbert I d 19741 1 19761) and by 
iLequeux et al. I (119791) . devoted respectively to the HII regions 
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in the Magellanic Clouds and in irregular and blue compact 
galaxies, derived a value of the helium to metals enrichment ra- 
tio of AY/AZ » 3 . From observations of extragalactic HII re- 
gions, iPagel et al.l (1 19921) fou nd that AY/AZ > 3, w ith a pre- 
ferred value close to 4, while iPeimbert et al. I (2000) obtained 
AY/AO = 3.5 + 0.9 from NGC346, the brightest HII region in 
SMC. A detailed study of a very large sample of spectroscopic 
ob servations of HII region s in blue compact galaxies performed 
bv lIzotov&Thuanl d2004l) led to AY/AZ = 2.8 + 0.5. However 
a reanalysis of the same data by Fukugita & Ka wasaki! ((2006) 
provides AY/AZ = 4-5, values that becomes as low as 1.1+ 1 
if stellar absorption is taken into account. 

From observations of galactic HII regions, such as the Orion 
nebul a, M8 and Ml 7, various authors found AY /AZ > 2 (see 
als o lEsteban at al.lfl 999: De harveng et al. 1 12000). As an exam- 
ple |Pd^e^nul20^lOTmd - Ay/AZ = 2.1 + 0.6 adopting 
a temperature fluctuation parameter t 2 = 0.037 ( recommended) 
and 3.8 + 1.1 adopting t 2 = 0.000. More recently. iBalserl (120061) 
derived AY/AZ - 1.41 + 0.62 from the analysis of the sole M17 
and S206, whereas taking into account also the HII regions be- 
longing to the Magellanic Clouds and metal-poor galaxies, she 
found AY/AZ = 1.6. The last author warns also on the system- 
atic underestimation of the measured h elium abundance cause d 
by a clumping in HII regions (see also Mathis & Wood 2005). 
Finally, a new estimate of the chemical composition of M17 
bv lCarigi& Peimbert (2008) yields AY/AZ = 1.97 + 0.41 and 
4.00 + 0.75 adopting t 2 = 0.036 and 0.000, respectively. 

A useful technique to empirically constrain the stellar 
helium to heavier elements enrichment ratio is the analysis of 
the chemical composition of planetary nebulae (PNe). These 
objects present the advantage of being numerous in a large 
range of metallicity and of allowing precise estimates of helium 
abundance. On the other hand, it is a well established result 
of stellar evolution that the chemical composition of low and 
intermediate mass stars, i.e. the PNe's progenitors, is strongly 
modified by the dredge-up episodes, which enrich the envelope 
with inner material previously processed by nuclear reactions, 
that is mainly with fresh helium and, in the case of III dredge-up 
in thermally pulsating AGB stars, also with heavier elements, 
mainly carbon. Thus, the chemical composition of PNe is not 
representative of the protostellar cloud and in order to derive 
the original interstellar abundances a correction must be applied 
which takes into account the evolution of the surface abundances 
of the progenitor star. In the early attempts to determine the 
AY/AZ ratio using the PNe in the late '70s, such a correction was 
negle cted and values in th e range 2.2-3.6 were obtained (see e.g. 



Table 3. Initial helium T„„ and metal Z,„, abundances provided 
by recent standard solar models by different authors (see foot- 
notes). M odels are c omputed adopting the heavy element mix- 
ture bv lAsplund et all (120051) . 



Dodorico, Peimbert & Sabbadin 1976; Peimbert & Serrano 



1980). In our knowledge, Chiappi ni & Maciel 1 11994) were the 
first to take into account the correction due to the evolution of 
the prog enitor star and t hey obtained 4 < AY/AZ < 6.3. More 
recently, iMaciel I (1200 ll) found 2.8 < AY/AZ < 3.6 adopting 
F P = 0.23 and 2.0 < AY/AZ < 2.8 with F P = 0.24. 

An important independent constraint on the stellar helium to 
heavier elements enrichment ratio is provided by the Sun, for 
which very detailed and accurate data are available. In order to 
derive the AY/AZ ratio from the Sun, the precise spectroscopic 
estimates of the solar photospheric chemical composition are not 
enough because they significantly differ from the original ones 
due to diffusion and gravitational settling. Thus, one has to rely 
on a standard solar model (SSM), that is, a stellar model of 1 
Mq which at the age of the Sun (i.e. ^4.56 Gyr) fits the solar ob- 
servables. As a result of such a procedure, the initial metallicity 
Zq and helium Yq abun dance are inferred. As l argely debated 
in recent years (see e.g iBahcall et al. I l2005bl: iBasu & Antia I 



Models 


Y 

L mi 


7 


AY/AZ 


Pisa" 


0.2532 


0.0137 


0.380 


B2005'' 


0.2614 


0.0140 


0.957 


G2006 r 


0.2570 


0.0135 


0.666 


S2009 rf 


0.2593 


0.0139 


0.813 


S2009b c 


0.2617 


0.0149 


0.919 



" iToenelli. P rada Moroni & Degl'Inno centil i2010h 
b Bahcall, Serenelli & Basu (2005) 
c iGuzi k. Watson & Cox] J2005D 

d ISerenelli et all ((2009) 

e Serenelli et al. (2009), mixture by Asplund et al. (2009) 



120081: IChristensen-Dalsgaard et al. 1 120091: ISerenelli et~ai1l2009l 
and references therein), the very good agreement between the 
SSM and the helioseismological constraints, mainly the sound 
speed profile, the extension of the convective envelope and the 
surface helium abundance, achieved at the end of the last century 
has been compromised by the new determinations of the metal 
abundances based on 3D photospheric models by Asplund et al. 
(200 5J), and only slig htly alleviated by the very recent release by 
Asplund et al. (2009). With these caveats, the initial helium and 
metal abundance of the Sun provided by the current SSMs allow 
to constrain the AY/AZ, once a primordial helium abundance Fp 
is chosen. Table [3] shows the results of SSMs computed by dif- 
ferent groups by assuming different solar heavy elements mix- 
tures. As can be easily seen, different authors provide values in 
good mutual a greement, values that t r anslat e in AY/AZ < 1, if 
F P = 0.248 bv llzotov. Thuan & Stasi I d2007l) is used. 

In principle the AY/AZ ratio can also be theoretically 
predicted by models of galactic chemical evolution. These 
predictions rely on several assumption on the initial mass 
function, star formation rate, star formation efficiency, stellar 
yields from intermediate and massive stars, ISM mixing; so 
even though they are quite powerful, due to the uncertain- 
ties on each of the aforementioned quantities they are still 
quite uncertain in predicting the expected helium to metals 
enrichment ratio and provide values which span a range 
of A Y/ AZ between about 1 and 4 ([Serrano & Peimbert 
1981; Timmes. Wooslev & Weaver! 1 1 995t iTosi I 1 1 996: 
Chiappini, Matteucci & Gratton "Tl997[ Fields & Olive 



199 8: iCarigil I2000t "^ Chiappini. Renda & Matteucci 2002; 
iRomano et al.l[2005t: ICarigi & Peimbertll2008l) . Among the pa- 
rameters that mostly affect the outcomes of chemical evolution 
models we mention some of the most important, such as the 
SN la rate, which strongly affects the iron enrichment, the star 
formation history and the initial mass function, since stars of 
different mass, which evolve on different time scales, contribute 
differently to the production of F and Z; the maximum stellar 
mass that enriches the interstellar medium is an other important 
parameter. Moreover, as far as the chemical enrichment of 
the interstellar medium is concerned, stellar models, which 
provide the yields to chemical models, are still significantly 
affected by poorly understood physical processes, such as con- 
vection, which affects the stellar yields through the efficiency 
of dredge-up phenomena and hot-bottom burning, mass-loss, 
which directly influence the enrichment, and rotation, which 
has a strong impact on the helium, carbon, nitrogen and oxygen 
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yields dMevnet & Maederl |2002). As a consequence, stellar 
yields provided by different authors still show not negligible 
discrepancies. Thus, the results of chemical evolution models 
depend on the adopted set of stellar yields ([F ields & Olive 
1 1998t ICarigi Il2000t IChiappini. Matteucci & Mevnetll2003l) . For 
example, it has long been known that the AY/AZ ratio grows by 
increasing the mass-loss rate of both intermediate and massive 
stars (see e.g. Ilben & Truran I [1978; R enzini & Voli I 119811: 

lPeimbertllT986l). 

Chiappini, Matteucci & Meynet (2003) predicted a value 
of the enrich ment ratio AY/AZ ~ 2.4, if the stellar mod- 
els by iMevnet & Maeder (2002), which take into account ro- 
tation, are adopted in their chemical evolution model, whereas 
a value of 1.5 is obtained, if the stellar computations by 
Ivan den Hoek & Groenewegen ( 1997), for the lo w and interme- 
diate mass stars, and IWooslev & Weaver! d 1995b . for the mas- 
sive ones, are used. The recent vers ion of the chemical evolution 
model bv lCarigi&Peimbertl d2008l) provides AY/AZ =1.70 and 
1 .62, if high or low wind yields are adopted, respectively. 



10. Conclusions 

The principal aim of this work was to test the reliability of the 
determination of AY/AZ by the comparison between low-MS 
stars and theoretical ZAMS models. A very fine grid of stellar 
models has been computed for many values of AY/AZ, [Fe/H] 
and masses adopting two different mixing-length parameters a, 
namely 1 .97 (our solar calibration) and 2.4, and two different sets 
of atmosphere models to transform luminosities and effective 
temperatures into magnitude s and color indices, the PHOENIX 
GAIA v2.6.1 models dBrott et alj|2005h and the ATLAS9 ones 
dCastelli & Kuruczll2003l) . respectively. 

A detailed analysis of the capabilities of the method and of 
the main uncertainty sources affecting the derived results has 
been performed by means of many numerical experiments on 
synthetic data set produced under controlled conditions and with 
precisely known properties. 

One of the main findings of the paper is that the inferred 
value of AY/AZ is quite sensitive to the age of the stellar sample, 
even in the case in which only very faint (i.e. My > 6 mag) MS 
stars are selected. By means of numerical experiment we showed 
that the lack of an age estimate of low mass field stars leads to 
an underestimate of the inferred AY/AZ of about 2 units. As a 
consequence, the face value of the helium-to-metals enrichment 
ratio provided by the recovery procedure applied to the solar 
neighborhood stellar sample must be corrected for this age-bias. 

Adopting our reference set of models (i.e. those with our "so- 
lar" calibrated a) transformed into the observational plane us- 
ing the PHOENIX GAIA v2.6. 1 model atmospheres, we found 
AY/AZ = 5.3 ± 1.4, 

Such a result has been checked against an independent and 
very accurate data set, that is, the low-MS of the Hyades cluster. 
This sample has the additional advantage to be unaffected by the 
age-bias, since its very young age (i.e. 500H-600 Myr) guarantees 
that the low-mass (< 0.9Mq) MS stars which constitutes the 
sample are essentially unevolved. The recovery method provided 
again AY/AZ = 4.75 ± 0.35, in perfect agreement with the result 
for field stars. 

To further check the consistency of this final result, we cal- 
culated isochrones with AY/AZ - 5 (the closest in our grid of 
models) and the measured [Fe/H] of Hyades. The good agree- 
ment between these isochrones and the Hyades MS, not only in 
the faint part belonging to the sample but also at higher luminosi- 



ties, is a further proof of the internal consistency of the recovery 
procedure. 

The effect of a change in the assumed efficiency of the su- 
peradiabatic convection (i.e. the a parameter) in the stellar mod- 
els used to build the ZAMS and in the adopted atmosphere 
models used to transform luminosities and effectives tempera- 
tures into magnitudes and color indices has been discussed, too. 
More in detail, the recovery method yield a nominal value of 
AY/AZ = 1.59 + 1.01 when adopting the set of theoretical mod- 
els computed with a = 2.4 and transformed into the observa- 
tional plane by means of PHOENIX atmosphere models and 
AY/AZ = 0.5 when adopting our standard set of models with 
a — 1 .97 but the ATLAS9 atmosphere models. These values be- 
come about 3.6 and 2.5, respectively, once corrected for the age 
bias. 

Our data have [Fe/H] determina t ions co min g from two dif- 
ferent sources, i.e. iNordstrom et alj d2004 and lTavlorl d2005h 
catalogs. The 4 stars we have in common from the two catalogs 
show a disagreement in the [Fe/H] scale betwee n the two cata- 
logs. This disagreement was already identified in lTavlorl {2005) 
and the author provides a Table to account for it and put the two 
catalogs on the same [Fe/H] scale. Nevertheless, even after the 
re-zeroing procedure suggested by iTavlorl d2005l) . we still find 
some disagreement in the final results on AY/AZ when the two 
catalogs are considered separately. Probably an even more ac- 
curate study of the zero points of metallicity determinations is 
needed. 
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